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Abstract 

We study a model of wealth dynamics [Bouchaud and Mezard 2000, Physica A 
282, 536] which mimics transactions among economic agents. The outcomes of the 
model are shown to depend strongly on the topological properties of the underlying 
transaction network. The extreme cases of a fully connected and a fully discon- 
nected network yield power-law and log-normal forms of the wealth distribution 
respectively. We perform numerical simulations in order to test the model on more 
complex network topologies. We show that the mixed form of most empirical dis- 
tributions (displaying a non-smooth transition from a log-normal to a power-law 
form) can be traced back to a heterogeneous topology with varying link density, 
which on the other hand is a recently observed property of real networks. 
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1 Introduction 



How networks self-organize into complex structures is one of the main inter- 
ests of modern statistical mechanics [1,2]. A large amount of data has been 
assembled which suggests that many different biological, social and techno- 
logical networks display non-trivial properties, strongly deviating from what 
expected in simple random graphs [2]. This has stimulated the proposal of 
theoretical models aimed at reproducing the observed topological features by 
means of basic mechanisms, invoking either growth processes [1,2] or 'hidden' 
underlying principles [3]. 

Understanding the topological properties of real-world networks is fundamen- 
tal also because they affect the outcomes of dynamical processes denned on 
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them: strong quantitative and qualitative effects of network topology have 
been highlighted in models describing synchronization [1], spread of epidemics 
[4] and wealth exchange [5,6,7]. In the present paper we focus on the latter 
topic by discussing how the form of the statistical distribution of the wealth 
of a set of economic agents depends on the topological properties of the trans- 
action network defined among them. Our analyses are performed using the 
model of wealth exchange proposed by Bouchaud and Mezard [5]. This prob- 
lem is interesting for various reasons: firstly, the outcomes of the model can 
be directly compared with empirical data [7,8,9,10,11]; secondly, wealth dis- 
tributions range in general between two typical forms which interestingly cor- 
respond to the two extreme cases of a fully connected and a fully disconnected 
network; finally, the problem is a first step towards the more general one of 
understanding how network topology affects the outcomes of a random process 
governed by stochastic differential equations. 



2 The Bouchaud-Mezard Model 



In the model of Bouchaud and Mezard [5] (BM in the following) the wealth 
Wi of an agent % {% — 1, ... AT) is assumed to be governed by the following set 
of stochastic differential equations: 

Wi(t) = rji(t)wi(t) + J2 JijWjit) - JjiWiif) (1) 



where the ?7j(t)'s are independent Gaussian variables of mean m and variance 
2a 2 (accounting for the intrinsic variation of wealth) and Jy = Jaij/N is the 
fraction of agent j's wealth flowing into agent i's wealth, which is directly 
related to the adjacency matrix of the underlying transaction network 
(which in the following we assume to be undirected: = — 1 if i and j 
exchange wealth and = a 3 -j = if they do not). Equation (1) is interpreted 
in the Stratonovich sense. 

The above model can be easily studied within the Fokker-Planck approach in 
the mean-field limit of a fully connected network (a^- = 1 Vi ^ j) [5]. In 
such the form of the equilibrium (long-term) wealth distribution p(x) 

(expressed as a function of the normalized wealth Xi = Wi/J2j w j) can be 
obtained analytically [5] and reads 



where f3 m f = 1 + J/cr 2 . The tail (large w limit) of this distribution clearly 
follows a power-law behaviour with exponent —1 — j3 m f. 
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In the opposite case of a network with no link (J^ = Vi, j), the sums ap- 
pearing in eq.(l) vanish and the process is purely multiplicative. The Central 
Limit Theorem then states that for large N the distribution of log{x) is nor- 
mal, or in other words the distribution of x is log-normal. A useful family of 
distributions is therefore given by 

expf-a/x) 

<* x i+p ( 3 ) 



which is almost undistinguishable from a log-normal one when (3 ~ and 
displays power-law tails of the form for (3 > 0. In particular, the case 

a = (3 — 1 = J I a 2 corresponds to the mean- field prediction (2). 

Real wealth distributions [8,10] always display power-law tails with exponent 
— 1 — (3 where (3 is usually near to, but often smaller than, one. By contrast, 
the mean- field limit of the BM model always predicts that f3 > 1. Another 
property of real data is in general a non-smooth transition from a left log- 
normal-like part of the distribution and a right power-law tail [7,8,9,10,11]. 
Both properties are not reproduced by the simple mean-field case, and in 
the following we ask if a suitable, and possibly more realistic, choice of the 
transaction network yields the desired outcomes. 



3 Random graphs, regular lattices and scale-free networks 



The first point of having an exponent (3 < 1 out of the BM model was 
already explored in [5], where the authors showed that when the network is 
a random graph [2] (characterized by a Poissonian degree distribution) (3 is 
systematically smaller than the mean-field prediction (3 m f. Clearly, when the 
connection probability p ~ all vertices are isolated and the log-normal-like 
distribution is recovered, while p — ► 1 yields the mean-field case. Therefore 
in the random graph case the distribution is well described by eq. (3), where 
(3 — > for p — > and (3 — > (3 m f for p — > 1 . 

We now show that a similar result holds when the network is a regular ring 
where each vertex is connected to its z/2 closest neighbours (that is, each 
vertex has degree z and the degree distribution equals P(k) = Sk z )- When z = 
2 one has a simple linear structure, and the log-normal-like (f3 ~ 0) behaviour 
is expected [5]. By contrast, when z — > iV — 1 the network approaches the fully 
connected topology and one expects (3 — > (3 m f. As we show in fig. la, this is 
indeed the case and the cumulative wealth distribution p>(x) = f£° p(x')dx' 
approaches the mean- field shape p> (x) oc x~^ m -t' as z increases (for simplicity, 
here and in the following we always set J /a 2 = 1 so that (3 m f = 2). Our results 
are in accordance with those of ref. [6], where the authors also introduce a 
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rewiring mechanism to obtain small-world networks [2] for various values of z 
in order to have a wealth distribution displaying a right power-law tail and a 
left log- normal part. They find that for a suitable choice of the parameters this 
indeed occurs, however the two parts of the distribution are always connected 
in a smooth fashion, in constrast with the non-smooth behaviour of real data. 




Fig. 1. Cumulative wealth distributions generated by the BM model on networks 
with N = 3000 vertices. The dashed line always corresponds to the slope predicted 
by the mean- field theory ((5=1 + J /a 2 = 2). a) Regular ring for various choices 
of the degree z. b) BA scale free network for various choices of the number m of 
vertices injected at each timestep. The solid line is drawn as a reference for the 
range of f3. 



We now turn to scale-free networks generated by the Barabasi-Albert (BA) 
mechanism [2], where at each timestep one introduces a new vertex and m new 
links connected to preexisting vertices chosen with probability proportional to 
their degree. The resulting degree distribution has the form P(k) oc k~ 3 [2]. 
We first generated a network with a given value of m and then let the BM 
model run on it. As shown in fig. lb, the cumulative distribution is always of 
the power-law type, with exponent (3 ranging from the quite small value 0.4 
to (3 = (3 mf = 2. 

The results discussed so far show that the wealth distribution generated by 
the BM model is strongly affected by the topology of the underlying network. 
The effect of the topology of the networks considered so far is however only 
quantitative, since essentially it consists in determining the exponent of the 
distribution (3). This allows to obtain values of the exponent smaller than the 
mean-field prediction, which was one of the necessary points we discussed at 
the beginning. However, in order to generate a mixed wealth distribution with 
the aforementioned non-smooth transition from a log-normal-like body to a 
power-law tail, deep qualitative effects of the topology are required, and more 
complicated networks must be considered. 
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4 Heterogeneously linked networks 



We now consider a very simple example to show how the mixed form can 
arise. Consider a network consisting of a fully connected sub-network with M 
vertices, plus a set of N — M completely isolated vertices. Clearly, in this case 
we expect the wealth distribution to be of the form p(x) = (M/N)p p i(x) + 
(1 — M/N)pi n (x), where p v \ and pi n indicate the power-law and log- normal 
predictions respectively. Then, by tuning the ratio M/N it is possible to vary 
the relative weights of the two contributions. The results are reported in fig. 
2a, showing that the mixed form is indeed observed. 

We can now introduce more complex models characterized by the same basic 
ingredient of a dense 'core' and a sparse 'periphery'. In the absence of empirical 
data on real transaction networks, all choices are quite arbitrary and here we 
simply report the results obtained for an octopus network having a random 
graph with M vertices as the core and a set of N — M 'tentacles' with single 
edges linked to randomly chosen vertices in the core. The results are shown in 
fig. 2b, and in a wide intermediate range the form of the distribution is indeed 
very similar to the observed one [7,8,9,10,11]. 
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Fig. 2. Cumulative wealth distributions generated by the BM model on a) mixed 
networks and b) octopus networks for various choices of the ratio M/N (dashed 
line: mean-field prediction j3 = 2). The number of vertices is N = 5000. 



5 Discussion 



We showed that a crucial ingredient in order to have a realistic wealth distri- 
bution out of the BM model is a heterogeneous link density in the network. 
Detecting regions with high link density is crucial in social networks, where 
clustered portions of the network usually indicate communities of tightly in- 
teracting people. More in general, the way the network can display heteroge- 
neously linked regions depends not simply on the degree distribution P(k), 
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but on higher-order properties such as the correlations between the degree of 
adjacent vertices [12] and the hierarchical dependence of the local clustering 
properties on the degree [13]. Indeed, in the present paper we showed that P(k) 
alone does not determine the desired properties of the wealth distribution. In 
order to relate the empirical form of wealth distributions to higher-order non- 
trivial topological properties would require the knowledge of the corresponding 
transaction networks. While this is a very hard task at the individual level, 
it is possible for the trade network of world countries, where the vertices are 
countries and links represent import/export trade relationships. Indeed, a re- 
cent study [14] of the topological properties of the World Trade Web (WTW) 
revealed the presence of degree correlations and of a hierarchical organization 
of the network. Interesting, these are exactly the properties that, according to 
our analysis, are expected to determine the peculiar shape of the corresponding 
distribution of the 'wealth' of vertices, which is this case is the Gross Domestic 
Product (GDP) of countries. Indeed, the GDP distribution of world countries 
closely resembles the mixed form shown in fig. 2b [7]. Future studies on the 
interplay between network topology of the WTW and the GDP dynamics, 
currently in preparation, will clarify this aspect in a more quantitative way. 
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